Marangoni shocks in unobstructed soap-film flows 
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It is widely thought that in steady, gravity-driven, unobstructed soap-film flows, the velocity in- 
creases monotonically downstream. Here we show experimentally that the velocity increases, peaks, 
drops abruptly, then lessens gradually downstream. We argue theoretically and verify experimen- 
tally that the abrupt drop in velocity corresponds to a Marangoni shock, a type of shock related to 
the elasticity of the film. Marangoni shocks induce locally intense turbulent fluctuations and may 
help elucidate the mechanisms that produce two-dimensional turbulence away from boundaries. 



Soap-film flows [lj have long been used to study two- 
dimensional (2-D) turbulence, a type of turbulence that 
differs from its three-dimensional (3-D) counterpart in 
crucial respects. For example, in 3-D turbulence the en- 
ergy may cascade only from larger to smaller lengthscales 
whereas in 2-D turbulence the energy may cascade in ei- 
ther direction [2| . Disparate directions of energy transfer 
result in disparate apportionings of the turbulent kinetic 
energy among the lengthscales of the flow Q ■ Besides the 
theoretical interest inherent in its distinctive characteris- 
tics, 2-D turbulence is relevant to the large-scale irregu- 
larities encountered in 2-D atmospheric flows — flows that 
are confined to two dimensions by geostrophic forces and 
a stratified atmosphere [3j. Examples of large-scale ir- 
regularities in 2-D atmospheric flows include hurricanes, 
typhoons, and the Great Red Spot of Jupiter [4|]. 

In the typical setup used to study soap-film flows [|| , a 
film hangs between two long, vertical, mutually parallel 
wires a few centimeters apart from one another. Driven 
by gravity, a steady vertical flow soon becomes estab- 
lished within the film (Fig. Q]). Then, the thickness h of 
the film is roughly uniform on any cross section of the flow 
and we write h = h{x), where x runs along the cen- 
tcrlinc of the flow (Fig. [T]). In a typical flow h w 10 /Ltm, 
much smaller than both the width w and the length L of 
the flow (Fig. [l}. As a result, the velocity of the flow lies 
on the plane of the film, and the flow is 2-D. Since the vis- 
cous stresses (and the attendant velocity gradients) are 
confined close to the wires, the mean (time-averaged) ve- 
locity u is roughly uniform on any cross section of the 
film and we write u = u(x). Thus, assuming incom- 
pressibility, h(x)u(x) equals the flux q per unit width of 
film and is independent of x for a steady flow. 

Analyses of steady flows have accounted for the grav- 
itational force, the inertial force, the drag force of the 
ambient air, and the drag force of the wires. Rutgers 
et al. Q have shown that the drag force of the wires is 
negligible as compared to the drag force of the ambient 
air and may be dropped from the equation of momen- 
tum balance. Then, a prediction can be made that in 
a steady flow the mean velocity is a monotonically in- 
creasing function of x and approaches a terminal veloc- 



FIG. 1: Typical setup used to study steady, gravity-driven, 
unobstructed soap-film flows. Axis x runs vertically along 
the centerline of the flow. Wires WL and WR are thin nylon- 
fishing lines kept taut by weight W. The film hangs from 
the wires; its width increases from to w over an expansion 
section of length I, then remains constant and equal to w over 
a measurement section of length L S> I- (The origin of x 
is at the top of the measurement section. L is the "length 
of the flow" and w the "width of the flow.") Reservoir RT 
contains a soapy solution which flows through valve V and into 
the film. After flowing through the film with mean velocity 
u(x), the soapy solution drains into reservoir RB and returns 
to reservoir RT via pump P. In our experiments, the soapy 
solution consists of ~ 2.5 % Dawn Nonultra in water; w = 2.5 
to 5.1cm; L — 1.05 to 1.39 m; and I — 23.5 cm. 



ity asymptotically downstream 0, @|. This prediction 
has not been tested, but it is thought to be in qualita- 
tive agreement with the few known experiments 
contrast to this prediction, in our experiments u(x) is a 
strongly non-monotonic function of x. 

To measure u, we use laser Doppler velocimetry 
(LDV). In Fig. [2] we show plots of u along the centerline 
of several representative flows. In each flow, u increases 
downstream up to a certain point whereupon it peaks, 
drops abruptly to a fraction of its peak value, then con- 
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tinues to lessen gradually downstream. 
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FIG. 2: Plots of the mean velocity u vs. x for three steady 
flows. The width w = 5.1cm for all flows; both the length 
L and the flux per unit width q change from flow to flow. 
L = 1.05m and q = 3.81CT 6 m 2 /s (A), L = 1.23m and q = 
4.910~ 6 m 2 /s (B), L = 1.39 m and q = 5.5 10~ 6 m 2 /s (C). 
Inset: Plots of u vs. x for three steady flows, w — 5.1 cm and 
q — 5.7 10 -6 m 2 /s for all flows; L changes from flow to flow. 

From the incompressibility condition (uh — q), the 
abrupt drop in u should be accompanied by an abrupt 
increase in h. To verify this abrupt increase in h, we 
light one face of a film with a sodium lamp and observe 
the interference fringes that form there. In Fig. we 
show a photograph of the interference fringes on the part 
of a film where u drops abruptly. The distance between 
successive fringes decreases rapidly in the downstream 
direction, signaling an abrupt increase in h. 

To verify the abrupt increase in h by means of an al- 
ternative technique, we put Flourescein dye in the soapy 
solution and focus incoherent blue light on a spot on 
the film. The spot becomes fluorescent, and we monitor 
the intensity of the fluorescence using a photodetector 
whose counting rate is proportional to h. In Fig. [5Jd we 
show plots of h along four cross sections of a flow. These 
cross sections lie on the part of the flow where u drops 
abruptly. The thickness trebles over a short distance of 
a few centimeters in the downstream direction. 

To explain our experimental results, we write the 
steady-state equation of momentum balance in the form 

phuu x = 2a x + pgh - 2r Q , (1) 

where p is the density, (-) x = d(-)/dx, a is the surface 
tension, g is the gravitational acceleration, and r a is the 
shear stress due to air friction. From left to right, the 
terms in ((TJ) represent the inertial force, the elastic force, 
the gravitational force, and the drag force of the ambi- 
ent air. Here we follow Rutgers et al. Q and use (as 
a rough approximation) r a = 0.3^/ p a p a u 3 / (x + I), the 
Blasius expression for the shear stress on a rigid plate 
that moves at a constant velocity u through air of den- 
sity p a = 1.2kg/m 3 and viscosity p a = 1.7 x 10~ 5 kg/ms. 
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FIG. 3: (a) Fringes over the part of a film where the velocity 
drops abruptly, (b) Plots of the thickness vs. z along four 
cross sections of a flow of width 2.5 cm and length 1.2 m. The 
cross sections are at x = 0.95 m (A), x — 1.04m (B), x — 
1.05m (C), and x = 1.07m (D). The large peaks near the 
lateral edges are due to backscattering from the wires. Fringes 
at a piercing upstream (c) and downstream (d) of the drop in 
velocity; fields of view = 5cmxl cm. 

To obtain an expression for a x , we argue that the con- 
centration of soap molecules in the bulk of the film re- 
mains constant in our experiments (because there is no 
time for diffusional exchange between the bulk and the 
faces of the film [1, Q). Then, the film is said to be in the 
Marangoni regime, and 2a x — —pU M h x [l(|, where Um 
is the Marangoni speed — a property of the film, indepen- 
dent of h, that quantifies the speed at which disturbances 
in h travel on the plane of the film 0, ■ By substitut- 
ing 2a x = —pUlfh x and h = q/u in ([T]), we obtain the 
governing equation 



In ^ we can distinguish two types of flow: a supercrit- 
ical flow in which u > Um and u x > 0, and a subcritical 
flow in which u < Um and u x < 0. We conjecture that in 
our experiments the flow is supercritical upstream of the 
drop in velocity and subcritical downstream. Consistent 
with this conjecture, for any fixed q the flow upstream of 
the drop in velocity remains invariant to changes in the 
length of the flow (inset of Fig. [2]) . 

To confirm that flows are supercritical upstream of the 
drop in velocity and subcritical downstream, we use pins 
to pierce a flow upstream and downstream of the drop in 
velocity (Figs. [3J: and d, respectively). Upstream of the 
drop in velocity the Mach angle w 50° (from Fig. [3J;) , and 
the local u = 1.83 m/s (from an LDV measurement); thus 
Um ~ sin 50° x 1.83m/s= 1.4m/s in our experiments. 
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Let us test the governing equation Q for one of our ex- 
periments. We adopt a value of U m and a value of q and 
perform two computations. First, we integrate @ down- 
stream from x = with boundary condition it(0) = u , 
where uq is the velocity measured at x = in the exper- 
iment [111 ]. This first computation gives a function u(x) 
that should fit the experiment upstream of the drop in 
velocity. Second, we integrate ([2|) upstream from x = L 
with boundary condition u(L) — ul, where ul is the ve- 
locity measured at x = L in the experiment [ill ] . This 
second computation gives a function u(x) that should fit 
the experiment downstream of the drop in velocity. We 
perform the same computations for each one of our ex- 
periments trying different values of q and Um-, an< i choose 
the optimal values of q and the optimal value of Um — 
that is to say, the values of q (one for each experiment) 
and the value of Um (the same for all experiments) that 
yield the best fits to the experiments (Fig. [4|. The op- 
timal value of Um, 1.48 m/s, is in remarkable agreement 
with our estimate from Fig. [3fc (1.4 m/s). The optimal 
values of q are in reasonable agreement with the experi- 
mental estimates for q [l2[ (caption to Fig. 2]). 
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FIG. 4: Plots of the computational u(x) (lines) and experi- 
mental u(x) (points) for ten different flows. The computations 
are for Um — 1.48 m/s and the values of q indicated below 
(the experimental estimates for q are indicated in parentheses 
[13]). w = 5.1cm for all flows, (a) Flows of length 1.05 m: 
(A) q = 5.710~ 6 m 2 /s (3.91(T 6 m 2 /s), (B) 2510 _6 m 2 /s 
(5.710" 6 m 2 /s); (b) flows of length 1.17 m: (A) 5.710" 6 m 2 /s 
(5.11(T 6 m 2 /s), (B) 7.41(T 6 m 2 /s (5.9 1(T 6 m 2 /s), (C) 
301(T 6 m 2 /s (7.510~ 6 m 2 /s); (c) flows of length 1.23m: 
(A) 6.1 1CT 8 m 2 /s (4.11(T 6 m7s), (B) 1410 _6 m 2 /s 
(5.310" 6 m 2 /s), (C) 311(T 6 m 2 /s (6.5 10" 6 m 2 /s); and 

(A) 1610" B m 2 /s (4.710- 6 m 2 /s), 
m 2 /s). 



(d) flows of length 1.39 m: 
(B) 25KT 6 m 2 /s (6.3 10~' 



We conclude that a drop in velocity signals a 



supercritical-to-subcritical transition and corresponds to 
a Marangoni shock. In theory the drop in velocity is in- 
finitely steep and may be said to take place at x = x*, 
where u attains the value of Um (and u x becomes singu- 
lar) in the subcritical flow (Fig. 2|) [13]. But in experi- 
ments the drop in velocity takes place over a finite span 
Ax whose magnitude appears to increase with q (Fig. [4]) 
and whose downstream edge is located at about x — x* , 
the theoretical position of the shock (a position which ap- 
pears to move downstream as q increases). Thus in our 
simple theory the shock is sharp whereas in experiments 
the shock is diffused over a finite span Ax. 

To understand the reason why our theory (which does 
not account for turbulence) cannot resolve the structure 
of the shock, recall that a shock must dissipate energy at 
a steady rate 14j . We argue (i) that the shock dissipates 
energy by powering locally intense turbulent fluctuations 
and (ii) that these fluctuations must extend roughl y oy er 
the same span Ax as the shock that powers them [15] . 

To test these arguments we use LDV to measure the 
root-mean-square velocity u rms along the centerline of 
a representative flow (Fig. [5]). From a comparison of 
Figs. [5^ and b, we confirm that the shock is accompanied 
over its entire span Ax by velocity fluctuations that are 
up to thrice as intense as the velocity fluctuations that 
prevail both upstream and dowstream of Ax. We con- 
jecture that intense velocity fluctuations can arise more 
readily where the mean velocity is higher; this may ex- 
plain why the locally intense velocity fluctuations — and 
the diffusive shock that powers them — are located on the 
supercritical side of the theoretical position of the shock. 

To verify that the velocity fluctuations are turbulent, 
we obtain the energy spectrum on the centerline of the 
flow for three cross sections: one upstream, one within, 
and one downstream of the shock (inset to Fig. [5Jd). 
(These cross sections are marked "A," "B," and "C" in 
Fig-El 3 -) The area under the spectrum is larger for cross 
section B than for cross sections A and C, confirming that 
the turbulence is more intense within the shock than else- 
where in the flow. Further, the slope of the spectrum at 
intermediate wavenumbers and the shape of the spectrum 
at low wavenumbers differ on either side of the shock, in- 
dicating that the spectrum undergoes structural changes 
as the flow traverses the shock. 

We have demonstrated the spontaneous occurrence of 
shocks in the soap-film flows that are customarily used 
in experimental work on two-dimensional turbulence. 
These shocks are dissipative and diffusive; they give rise 
to fluctuations independently from the boundaries, with 
a strong but circumscribed effect on the spatial distribu- 
tion of turbulent intensity; and they alter the structure of 
the turbulent spectrum downstream from the shock. We 
conclude that the presence of shocks should be factored in 
in the interpretation of experimental measurements, and 
submit that shocks may furnish a convenient setting to 
study localized turbulence production in two dimensions. 
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FIG. 5: (a) Plots of the computational u(x) (lines) and the 
experimental u(x) (points) for a representative flow. Aa; is the 
span of the shock, (b) Plot of the experimental uj m8 (an index 
of the energetic contents of the velocity fluctuations) vs. x for 
the same flow. Inset: energy spectra at the centerline of the 
same flow for the cross section at x = 0.60m (A), x = 1.17m 
(B), and x — 1.30m (C). These are the cross sections marked 
"A," "B," and "C" in part b. The spectra are log-log plots of 
the energy density E (m 3 /s 2 ) vs. the wavenumber k (1/m). 
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